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We consider finite-size effects on the radiative energy loss of a fast parton moving in a finite 
temperature strongly interacting medium, using the light cone path integral formalism put forward 
by Zakharov. We present a convenient reformulation of the problem which makes possible its exact 
numerical analysis. This is done by introducing the concept of a radiation rate in the presence of 
finite-size effects. This effectively extends the finite-temperature approach of AMY (Arnold, Moore, 
and Yaffe) to include interference between vacuum and medium radiation. We compare results with 
those obtained in the regime considered by AMY, with those obtained at leading order in an opacity 

■ expansion, and with those obtained deep in the LPM (Landau-Pomeranchuk-Migdal) regime. 
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Q '■ I. INTRODUCTION 

: 

The physics of strongly interacting matter in extreme conditions of temperature and density continues to be one of 
the most active and fascinating fields of contemporary subatomic physics. It straddles nuclear and particle physics, 
O ,' and enjoys an immediate connection with the science of dense stellar objects. The liveliness of this discipline owes 
in no small part to the vibrant experimental programs and active theoretical support that characterize it. In this 
context, it is fair to write that the Relativistic Heavy Ion Collider (RHIC) continues to be a center of activity, and that 
this facility now boasts several compelling results [lj]. For example, the measurement of a strong suppression of hard 
!— i. partons in nucleus- nucleus collisions has been spectacular Q, confirming - and even exceeding - early expectations 
[3|. In addition, the empirical success of the hydrodynamical modelling of soft hadronic observables has been an 
important factor in popularizing the concept of a strongly coupled quark-gluon plasma (sQGP) [J]. 

■ The experimental measurements of the quenching of hard QCD jets have stimulated the theory community and 
\ several approaches based on perturbative QCD (pQCD) have been put forward to explain them. The formalisms 

differ in the details by which they treat the QCD medium, but all include the coherence built in multiple scatterings: 
The so-called Landau-Pomeranchuk-Migdal (LPM) effect @. Recently, some work has been devoted to comparing 

\q . the different theoretical approaches, to highlight the similarities, and to point out the differences [6j. Importantly, the 
medium in a relativistic heavy ion collision is a dynamical one: It rapidly evolves with time. It is therefore important to 
separate the consequences of the dynamical evolution scenario from those related to basic and fundamental differences 
^ ■ in the theory [7] . In order to highlight as clearly as possible the different assumptions and ingredients going into the 
different theory approaches, it has been useful to imagine an idealized setting where the basic formalisms can simply 

• '-J i be compared without the complications of realistic dynamics: The "QCD brick" . The brick is a slab of strongly 
interacting partonic matter, in thermal equilibrium at a temperature T. The goal then consists of analyzing the 
energy loss of a fast parton, propagating in, and interacting with, the medium. With this philosophy in mind, we 
concentrate in this paper on AMY, a formalism which stems from finite-temperature field theory and treats the QCD 
medium dynamically. Its basic premises were established in a series of papers written in the last decade [1, Q , and 
some of its phenomenological consequences have also been investigated 0, [1(1 [H| ■ The original formulation of AMY 
was in momentum space, and effects owing to the finite size of the emitting region were therefore a challenge to include. 
This paper follows the techniques of the light-cone path integral formalism elaborated by Zakharov [12. Il3j and which 
is known to lead to results similar to those of Baier-Dokshitzer-Mueller-Peigne-Schiff (BDMPS) [lj] and incorporate 
these effects into the AMY description. The effects of a finite formation time on inelastic radiation rates and on the 
time-dependent parton momentum distributions are explored, in order to eventually assess the consequences of this 
additional consideration on the phenomenology of RHIC and on that of the LHC. 

Our paper is organized as follows: In the next section, the basic theoretical building blocks of our approach are laid 
down. In what regards radiation rates, we show that in well-defined limits our results coincide with those obtained in 
other treatments. Next we explore numerically the effect of the finite formation time on radiation rates, and then on 
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the parton distribution functions. We explore the dynamics of the radiated gluons, and then summarize and conclude. 



II. FORMALISM 
A. Starting point 

The starting point here will be Eq. (4) of Ref. [l3| expressing the total probability for parton a, produced at time 
t = with energy p, to emit a pair b, c of bremsstrahlung partons with energies fc, p—k: 



dP b a c P^\x) x 
dk irp 



f°° f°° d d 
Re / dt 1 / dt 2 ———[K(t 2 ,*n;t u y 

Jo Jti dx d y 



) - (vac)] x=v=0 . (1) 



For the benefit of the reader, this equation will be re-derived in the next section. In Eq. (JXJ) , K(t 2 , x; t\, y) is 
\p/(2k(p— k))] 2 times the propagator associated with the light-cone Hamiltonian 

H = 8E(p)-iC 3 , (2a) 

6E(P) = 2k{Pk) + 2fc + V (2t) 

»^ / \ Cf,+C c — C a . . C a +C c — Cb . fc C a -\-Cb—C c f P~k . , . 

c 3(x) = w 2 (x ) + u 2 (-x) + v 2 ( x, (2c) 

2 2 p 2 p 

which acts on wave-function ?/>(x) on the transverse plane; m a ^, c and C a ,b,c are the respective masses (including 
thermal corrections) and group theory Casimirs of the particles involved (Ca — 3 for gluons and Cp = | for quarks in 

QCD); Cfi>2(x) is the dipole propagation amplitude per unit length for a qq pair. P^°\x) are the DGLAP splitting 
kernels (l5| for the process a — > be, with x = - the longitudinal momentum fraction of particle b: 

5 2 C F i±^, q^gq 
P f\ x) = I 9 Ca } ' , 5 -> 55 (3) 



2«? 2 jV F T F [a; 2 + (1-z) 2 ] , g ^ qq 
2 l+(l-xY 

* x e->je 

Final states be and cb are equivalent and should be counted only once. In particular, only the region x < \ is to be 
included in g -> g.g. For QCD, C F = f , Ca = 3 and N F T F = §. 

At this stage, it is appropriate to observe that the formula contains two time integrations, which, ultimately, 
correspond to two emission vertices - one in the matrix element and one in its complex conjugate. Since no emission 
occurs before the jet is produced, both integrations begin at t = 0. But since the radiation can be emitted at any 
time afterwards, they extend to infinity. Hence (JT|) contains two non-compact time integrations, which, in general, 
have to be performed numerically. 



B. Re-organization 



We re-arrange Eq. ([T]) in two steps in order to make it more suitable for our analysis. Mainly, we wish to remove 
e non-compact time integrations. First, we Fourier-transform it to transverse momentum space, with f k = J ^rp, 

dP b a c P b T\x) 



dk up 

In Fourier space, C3 becomes the Boltzmann-like collision operator 



Re / dti / dt 2 / q-p[Jf(£ 2 ,q;£i,p) - (vac)] 

JO Jti Jq.p 



cV(p) = l c ^{ Cb+C 2 C > (p)-V>(p-q)] 



+ c a +c a mp) _^k }] 
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and C(q) is related to the dipole amplitude via v 2 (x) = J q C(q) (l — e jq x ). This function can also be described as 
the Casimir-stripped elastic collision rate for a hard particle: 

c(q) = Q (27r) 

No confusion should arise between the function C(q) and the Casimirs C s . 

The second step is to make the time integrations compact. At large times, there are no collisions and K depends 
on time only through phases, which can be integrated exactly. We use this when integrating by parts: 



/•oo />oo 

/ dt 2 K(t 2 ,q]t 1 ,p) = / dt 2 
Ju Jt, 



d / e - i5E (i) t 2 



dt 2 V -iSE(q) 



««(«0*»A-( t2jq ;t 1)P ) 



-if (oo, q; t lf p) 1 , K(t 1 ,q; t u p) 



-i6E(q) v '™ -i8E(q)' 



OC 



+ J t dtj^C(t)K(t,q;t u p). (5) 

The first, rapidly oscillating, term on the second line disappears, owing to the conventional convergence factor e~ et 
implicit in Eq. (pj. The second term is closely related to the vacuum (DGLAP) radiation, and for the moment we 
simply cancel it against the vacuum subtraction. In fact, this cancellation is not exact and a remainder associated 
with the thermal masses will remain; we postpone its discussion to section HVl For now, only the last term survives, 
giving 



xRe / dt f dh f ■£^C(t)K(t,q.;t 1 ,p). 
Jo Jo Jn,r, oE(q) 



dk up Jq Jo 7 q! p 5E(q) 

In this form, it is natural to take a derivative with respect to t to define a radiation rate: 

druv = p b T(*) x Re r dti r c(W , q;tl , p) . (6) 



dk irp ' J Jq p SE(q) 

This is the formula which we shall study in this paper. It defines a radiation rate in terms of a single, compact, time 
integration that may conveniently be computed numerically. By construction, its content is fully equivalent to that 
in the BDMPS-Z probability, Eq. (JXJ) ; its time integral reproduces the latter. Details of its numerical implementation 
are given in Appendix [X] 

It is important to mention that a radiation rate is never uniquely defined. Indeed, the emission time of a given 
quanta is itself ambiguous up to a so-called formation time. Different prescriptions may lead to different time- 
dependent rates all associated with the same time-integrated probability. The present prescription is to locate the 
emission at the time t, which corresponds physically to the time of the last collision that influences the emission. 
This is a reasonable choice, as it is located between the two emission vertices. It is also a particularly natural and 
convenient choice because it causes the medium-induced radiation to stop outside the medium, that is, T is zero when 
C vanishes. In particular, for the brick problem, T is zero after exiting the brick. Nevertheless, all possible interference 
effects (including LPM interference), and memory effects, are included through the non-local kernel K. 

III. PHYSICAL INTERPRETATION OF THE FORMULA 

We now provide a heuristic derivation of our formula ([5]), as promised. This will help highlight its underlying 
assumptions and regime of validity. In this section we also discuss various analytical limits which are related to 
approximations used in the literature. 

A. Heuristic derivation of Eq. (flj) 

Numerous derivations of formulae similar to (JXJ) exist, exploiting one formalism or another 0, [l^, [3 [3 E3- We 
believe that the assumptions to be presented here, which did not previously appear with such generality, suffice for 
its validity. 

The four assumptions we will need to make are: 
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FIG. 1: A typical matrix element for gluon emission along with its complex conjugate. The initial blob represents the hard 
process that creates the hard parton a, drawn as a solid line, while the other blobs represent an arbitrary number of plasma 
particles. There is a monotonic time flow from the left to the right of the diagram which is a consequence of the eikonal 
approximation. The collision sub-diagrams, or "rungs", are assumed to have negligible horizontal extent. Only region B 
actually contributes to the emission process. 



1 The three "hard" particles a, b and c propagate eikonally. 

2 The opening angle of the radiation is small and its emission vertices are described by the leading-order DGLAP 
expressions. 

3 The elastic collisions against the plasma constituents are instantaneous compared to the formation time of the 
radiation. 

4 The transverse size of the radiation, at its formation, is small compared to medium scales. 

The meaning of assumption 1 is that the change in the propagation direction of any of the hard particles, over a 
formation time, is small. This does not necessarily mean that the opening angle should be small, only that it should 
stay constant with time. This will always be true provided E/T is large for all three particles, and corrections will 
come in powers of T/k, k denoting the smallest of the energies of the particles. 

The meaning of assumption 2 is that the emission itself is mediated by a leading-order DGLAP vertex. This does 
require a small opening angle, although a parametrically small angle may not be needed. This also requires that loop 
corrections be small. The relevant loops are presumably controlled by a s at the scale p± which enters in (|6]). An 
argument for this can be found in Ref. [HI, section VI B. The values of p± that appear in RHIC and LHC contexts 
will be discussed below. Incorporating these loop corrections would be difficult at present. 

The meaning of assumption 3 is that the coherence time of the force exerted on the hard particles is short compared 
to the time scale of the emission process. This coherence time can be estimated as ~ 1/VkT for the hardest elastic 
collisions up to ~ 1/T for plasma-scale ones (~ 1/gT for a typical soft collision if the plasma is weakly coupled). On 
the other hand, formation times scale like 

k k 1 / 2 
t ~ Bethe-Heitler regime, t ~ ^3/2 ' LPM regime. 

Formation times increase with a power of energy, while coherence times do not. Hence this assumption will always 
be valid up to T/k corrections. 

Assumption 4 means that the medium can be treated as homogeneous in the transverse plane. Note that no 
homogeneity assumption is made in the longitudinal direction. The case of a general, time-dependent medium is dealt 
with in Appendix |A"1 

Within these assumptions, the most general Feynman diagram contributing to a squared matrix element for 
bremsstrahlung is of the "ladder" form shown in Figure Q] All such diagrams contribute at the same order and 
have to be summed. The eikonal approximation implies a monotonic time flow along the horizontal hard lines a, b 
and c, and only two hard vertices connect those lines because of assumption 2. (One is in the matrix element, the other 
is in its complex conjugate.) The vertical rungs describe real and virtual elastic collisions that occur against plasma 
constituents. We allow the rungs to be arbitrarily complicated, multi-loop, objects - we will only need to assume that 
they have no horizontal extent. Because of this, all diagrams are cleanly separated by the emission vertices into three 
regions A, B and C. As we shall now review, only region B contributes to radiation probabilities. 

Consider, first, elastic collisions in region C ("elastic" refers only to the particles in the jet; the plasma particles 
could be scattered very inelastically). These collisions will change the momenta of the outgoing particles and typically 
reduce their energies by a relatively small amount, which will modify their momenta measured in the detector. To 
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discuss the radiation probability, the proper thing to consider is the integrated measurement over the hnal momenta. 
This is represented symbolically by the operator 

/ |P6jPc> ® <P6,Pc| = / d 2 x b rf 2 x c |x fc ,x c )(x b ,x c | (7) 

-'Pb.Pc J 

being convoluted on the right-hand side of the diagram in Figure [TJ In the eikonal approximation, this operator is 
time-independent. To put it as it is, "elastic collisions are elastic". Hence the sum over states can be performed at 
time t2, making region C irrelevant for the radiation probability. For the same reason, elastic collisions occurring in 
region A will not affect the emission probability, as they can be absorbed by a change in the jet axis. 

Thus only region B contributes. In the eikonal approximation, the object that is evolved in that region is a 
wavefunction of the form |qt,q c ) ® (<la|, depending on three transverse momenta. Equivalently, it depends on three 
transverse positions. This describes the eikonal amplitude times its complex conjugate in region B. Elastic collisions 
modify all three momenta but preserve the combination (q;,+q c — q a ). In fact, this combination vanishes at all times. 
Furthermore, it is possible to "gauge" to zero q a (or q b or q c ) by a rotation of the beam axis Q. Provided the opening 
angles are small, this rotation reduces to linear shifts Jq, ~ pf , which commute with both the eikonal expression for 
the energy ~ ^ qf /p, and with the collision term. Thus this dynamical gauge fixing does not introduce any new 
term, and the wavefunction really depends on a single transverse momentum. The evolution of this wavefunction will 
now be expressed in terms of the light-cone Schrodinger equation @ . 

Choosing to gauge q a to zero so that q = q b = — q c , the light-cone energy of the state \qt>, q c ) <g> (q Q | is 



Eh+Er — E„. 



2 2 

pq m t 



2k(p-k) 2k 2(p-k) 2p 



This gives SE in the Schrodinger equation @. 

The second ingredient in that equation is the collision operator. Quite generally, within the eikonal and instanta- 
neous approximations, the vertical rungs in Figure ([1]) reduce to a potential 

— iv(x, — x) 
P 

added to the Schrodinger equation. This is imaginary, corresponding to a collision rate. Here we have used assumption 
4, homogeneity in the transverse plane, to reduce it to a function of two separations, but we could still be left with 
a quite complicated function of two arguments. To proceed further, we need to make the extra assumption that v is 
well-described by the dipole-factorized form (p?c]) : 



C b +C c -C a C a +C c -C b .k . C a +C b -C c .p-k 
«a(x) + «a(-x) + v 2 (— x), 

for some function i>2(x). We believe that this is reasonable. In particular, it was shown in pl| that 0(g s ) effects at 
weak coupling do preserve this form. Perhaps, then, the dominant higher-order effects also do. 1 In momentum space, 
this assumption implies that the collision operator takes the form Q for some function C(q). 

Since little or no experimental data presently constrains C(q), and since its theoretical determination is so uncertain, 
at this stage it is probably best to think of the whole curve C(q) as a phenomenological parameter. We shall not have 
much more to say regarding its choice in this paper. 

The final ingredients which enter the diagrams are the hard vertices, which give the DGLAP factor in (fTJ. 

This forms the physical basis of the BDMPS-Z formula. Diagrammatically, to obtain our equation ([6]), one introduces 
a new time variable t which corresponds to the time of the last scattering before . The integral over t<i can then be 
performed analytically, since there are no collisions between t and i2, and time t becomes the emission time in (|6]). 



B. Leading order in the opacity expansion 

Gyulassy, Levai and Vitev [l6l.l2l| perform an expansion in the number of in- medium collisions, N, occurring during 
a formation time. Recently, it was shown that this expansion is equivalent to the series expansion of the BDMPS-Z 



Incidentally, color structures of the dipole type have been observed to survive NLO corrections in the seemingly unrelated context of 
soft wide-angle radiation |2(J . 
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formula (p} in powers of C [22|]. Furthermore, its numerical evaluation appears to converge up to N = 9. In that 
sense the GLV expansion is equivalent to the BDMPS-Z formula. Our equation ([5]) provides an exact and numerically 
tractable implementation of BDMPS-Z, and its series approximation is thus guaranteed to agree with the GLV series. 
For what follows it will be useful to highlight the N = 1 limit, which will be shown to work very well in certain 
regimes. 

In the N = 1 approximation, K reduces to the free Schrodinger propagator, reducing (j6|) to 



dk 



4nk 2 (p-k) 2 
4nk 2 (p-k) 2 



SE(p) 



pRci j dt\e 
o 



qp 



1 - cos(<5£(p)£) 

c(q) sm x 



-i<5 E (p)(*-*i) 



p 2 p (p q) 

5E(p) 5E(p-q) 



2 terms 



(8) 



The "2 terms" in the square bracket refer to the last two terms in If we drop the thermal masses in 5E and 
make the soft gluon approximation (small k), under which the bracket simplifies considerably, this precisely reduces 
toEq.(113) of Ref. [Hf. 2 

In the following, we will refer to ([5} (including thermal masses and without making the soft gluon approximation) 
as the N = 1 approximation. It is straightforward to evaluate numerically, and it will be seen to be valid at small 
times. 

A word about the collision kernel is in order. GLV conventionally use the "static scattering centers" approximation 
in which 



C(q) 



9s E CjUj/di 
(q 2 +Ai 2 ) 2 



(9) 



with rii, Ci and di being the number density, Casimir, and the representation dimension of various plasma constituents. 
This is a static Yukawa potential, with screening scale fj,. In general, this gets modified by recoil effects, which enhance 
the elastic cross-sections. As has been verified recently in great detail [23| (in the N = 1 case, but with a plausible 
argument for the general case), a simple change of C(q) to Eg. (fT2"|) given below will fully account for this effect. 
In fact, this particular finding of [23| is a consequence of the general argument we have presented above, since the 
assumptions made in (23j are not more general than the ones we have considered. This result is perfectly in line with 
the viewpoint advocated in the present paper, which is to view C(q) as a phenomenological input to the formalism. 
In this work we will only use the collision term (p~2]) below. 



C. AMY 



Some work has been done, relating the finite-temperature approach of Arnold-Moore- Yaffe (AMY) to that outlined 
in section III A[ in Ref. [I?} and independently in Ref. [24] . Other recent work on the AMY formalism includes [25| . 
The computation in the present subsection is very similar to those in Refs. [13, HH . In both references, the key step 
was to re-express BDMPS-Z in terms of a rate. This is also the idea behind ©. It is tantamount to formulating 
AMY in coordinate space. AMY consider bremsstrahlung from an homogeneous medium of infinite length and neglect 
effects associated with the jet production time t = as well as effects resulting from changes of the medium during a 
formation time. For the brick problem, this means that AMY should be valid in the large t limit. 

Our formulation connects very naturally with that of AMY. In the large t limit © becomes 

If-iS^^//^ ! ^ iJ ^L^lA)° c ° K ^ h -* y (10) 



2 GLV write all results in terms of a scattering length A and a q-distribution of collisions |U(q)| 2 normalized to unity. However, only the 
combination (27r) 2 | : u(q)| 2 /(AC8) appears, which is what we call C(q). With this substitution (and Ldz — > dt), their Eq.(113) literally 
becomes 



J pa 



C a C g g s 



k 

1? 



,(k-q 2 1 
-cos(^ —t) 



dkd 2 k 2n 3 k J 7 q v ' (k-q) 2 L (k-q) 2 
Upon stripping off the dt integration to define a rate and integrating over k, this reproduces ((Sj with p 
approximations are made into the latter. 



k q, provided the said 
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Now we observe that 

f( P ) + JL = f( P ) = 4fc2 ^ fc ) 2 f ^ f qi r(o,q;* 1)P ) 

has the same real part as f (p). By time-translation invariance, f(p) obeys the steady-state equation 

[iSE + C]~f(p)=p (11) 

which suffices to determine it. Eqs. (ITU)) and (fTTT) are precisely equations (fTP)) and (fTTj) appearing in [Pj. 3 Hence the 
rate which we have defined agrees precisely with that defined by AMY when t £f orm . 

It is sometimes stated that the AMY treatment, being rigorously justified in thermal field theory in the limit of 
small coupling, can be valid only for weakly coupled plasmas. This is incorrect. We see that this treatment will always 
be valid at times larger than a formation time, within the rather general assumptions discussed in section UlI Al On 
the other hand, a rigorous application of AMY does suggest the use of 

c(q) = J^llir with ml = ^ T2 in QCD ' (12) 

as the collision kernel. This expression is obtained at leading-order in the coupling in thermal field theory. (This was 
first obtained in this form in [26j.) This is the only place where a weakly coupled plasma is actually needed, when 
k T and the conditions of section IIII Al are satisfied. Otherwise, if the plasma is not truly weakly coupled but the 
jet still is, corrections to (fT2|) are to be expected. 



D. The deep LPM regime: The harmonic oscillator 



At asymptotically high energies p,k 3> T the deep LPM regime is realized, in which a large number of collisions 
occurs during a formation time. This motivates a diffusion approximation in which the Hamiltonian becomes 

it - PP 2 , 9 d 

H = m ^ + 2k(^k) + -d^W (13) 



where ih e s = + 2(p-k) ~ ~5p m ^ nc effective diffusion coefficient q 3 is related to the quark momentum broadening 
coefficient through q 3 = ^ c b +c.-c a + c a +C c -C b + c a +C b -C a (p^kf_ ^ where ^ . g defined as 



This Hamiltonian is that of a harmonic oscillator with complex frequency uiq — y 2 k(p~ P k) ■ This means that it can be 
solved analytically. For a homogeneous brick the exact solution is 

K(t,ci;t u p) = —3- TToTp- ; 7— — -exp 1 — : — — im e z{t-t x ) 

4fc J (p— ky zttiluq sin wo fa) \ 2kyp— fcjwo smwo(t- si J 

Upon performing the p integration in the rate ([6]) and changing variable to u cx q 2 , we obtain 
d^tc , n Q c (0) (^) x p ~Ht)m eS f 1 ^ e-^' sin w f f 00 udu e~" 



dk 4n 2 k(p-k) w Jo cos 3 w t' J a (u + i^ sinuiot ' ) 3 



It can be shown that denominators are always non-singular. 

The brick result, Eq. (|15p . will suffice for the present paper. Possibly, similar formulae for evolving media could be 
derived following the methods of [24]]. We will refer to (|T5|) as the "multiple soft scattering" approximation. 



3 The precise translation is that h in corresponds to pp here while F corresponds to 2p f . 
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It is unusual to express results obtained within this approximation in terms of a rate, and, for this reason, t|l 5|) may 
appear unfamiliar to many readers. We insist to use a rate in order to facilitate its comparison with ([5]) later on. To 
see that (fT5|) is indeed totally equivalent to the usual results, consider, for instance, its limit as m e g — > 0: 



dV a 
aL be 

dk 



8ir 2 k(p-k) 



q 3 (t) Re / dt' 



cos 2 (u;ni) 4ir 2 p 



x Re [—u>o tan u>ot] . 



(16) 



This is the time derivative of the familiar BDMPS formula = rh l^i^' In | costal, hence it is completely equivalent 
to it. More generally, when masses are included, we expect complete agreement with the formulae of [l2j . 

Physically, this approximation is justified when a large number of soft collisions occur during a formation time. 
While this is formally the case in the high-energy limit, the corresponding expansion parameter is only l/log(£') 
HP . 1 We will nevertheless verify that this approximation is accurate in the regime of high energy and sufficiently 
large time. At times smaller than a formation time, however, this approximation is not well motivated since a large 
number of collisions may not be reached. This was emphasized in |27|, H^|. In particular, in (2?J it is explained that 
formation time effects are only well described in a relatively restricted time window of parametric size y/log(E /T). 
This quantity is not large at RHIC (nor will it be at the LHC), and its smallness will become apparent numerically 
below. 



o(0) 



(x) 



IV. MULTIPLE EMISSIONS 

So far we have discussed how to calculate the rate at which a single quanta is radiated. One advantage of a rate 
formulation is that it is obvious how to account for multiple emissions: One simply writes down a "rate equation." 

In this section we keep the notation general, in anticipation of situations in which the local plasma frame will be 
evolving with time. Hence we work entirely in the lab frame. In general, the evolution of the distribution of quarks 
and gluons within a jet may be described within the Boltzmann framework 



dn q {p) 
dt 



dn 9 {p) 
dt 



= / dk 



n q (p+k)- 



ir ' ' ' dk 
[elastic] + [vacuum], 

dk 



dk 



-0(p-k) 



dT qg (p+k,k) dTS q (p+k,k) dTi qg (p,k) 

^Tk + 2n ° (p+k) ^ L Tk n «W— a 

.urn] , 

n q (p+k) — ^— Vn g {p+k)- 1 <-»' ■ '"■ 



-[elastic] + [vacuum], 



dk 



, dT g „Jp,k) dT 9 nn (p,k) 
~ ng{p) I dk e{p ~ k) + 0{p-2k) 



dk 



(17) 



where the k integrations run from to oo. The rate dT% c /dk(p, k) for emission of parton b with energy k from parent 
parton a with energy p is as defined by ([6]). In an evolving medium its computation is described in the Appendix 
IA1 Our notation follows that of [H, although the integration region is different; @ and other AMY-based treatments 
typically include absorption processes from the medium. The reason we omit these processes will be explained below. 

The "elastic" component accounts for the energy lost by collisions against plasma constituents. It has been exten- 
sively discussed in the literature (29j j . and we will not do so here. We only mention that this component is important 
when discussing observables such as Raa [30]. 

It is obvious that a vacuum component has to be included. This component, however, deserves further comment. 
We must recall what is already included in dT^ c /dk(p,k). This is the vacuum-subtracted rate © which includes all 
interference effects between vacuum and medium radiation. What is left is precisely the amount subtracted in ([5]) 
corresponding to the vacuum radiation, which may be written 



1 7->a,vac 
dP bc 

dk 



2irp 



2k(p-k) J 



dte-*fo SE 



(18) 



In the absence of thermal and vacuum masses, the integral is J q ^ which is the standard collinear logarithm for 
vacuum radiation. But because thermal masses are nonzero, (| 18(1 describes the probability for vacuum radiation with 



4 The expansion is in logarithms as opposed to powers essentially because collisions which are not soft always contribute. This is reflected 
in the logarithmic dependence of the parameter q Q14t on an ultraviolet cutoff; collisions near the cutoff are, by definition, not soft. 



Technically, the cutoff behaves like q n 



(QE) 



1/4 



18 
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FIG. 2: (Color online) Radiation rate for a 3 GeV gluon from a 16 GeV parent quark as a function of travelled length since 
the birth of the jet. The medium is a uniform brick of QGP at T = 0.2 GeV (panel (a)) and T = 0.4 GeV (panel (b)), with 
a s = 0.3. 



a slight suppression due to the thermal masses. This suppression is known as the (thermal) Ter-Mikaelian effect and 
it has been analyzed in the RHIC context in JUH^]. We will not discuss it further here. 

The single-emission probability (|18p is to be exponentiated so as to maintain a reasonable ordering of events; by 
our notation we do not wish to suggest that the vacuum radiation can be naturally handled as a rate. Some vacuum 
radiation will occur before or overlap with the earliest medium radiation, and some will be fragmentation radiation 
occurring afterwards in the confined phase. We hope to return to this question of ordering in a future work. For 
the application we have in mind in this work, we return to our idealized thermal medium, and concentrate on the 
radiative component dT^ c /dk. 
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FIG. 3: (Color online) Same as Figure [2] but for a 8 GeV gluon radiated from a 16 GeV quark. 



V. RESULTS AND DISCUSSION 



In this section, we calculate the scattering rates in different approaches and compare the results. Specifically, 
we compare results obtained with the formalism presented here, in which the formation time/length is included 
explicitly, with those obtained to first order in the opacity expansion, with AMY, and in the multiple soft scattering 
approximation. The first such comparison is shown in Figure [2j 

AMY is seen to be valid for large times, and our results do tend to that limit as r — > oo: A satisfying consistency 
check. There is a slight overshoot at a finite time, followed by asymptotic convergence, which we interpret as a gradual 
setting in of the LPM suppression. The left panel of Figure [2] shows the differential rate at which a 3 GeV gluon would 
be radiated off a 16 GeV quark in a "brick" of equilibrated quark-gluon plasma at a fixed temperature of T = 0.2 
GeV. The right panel represents results for similar requirements, but for a temperature of 0.4 GeV. For the coupling 
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FIG. 4: (Color online) "Instantaneous" radiation spectra k 2 dT/dk in a T = 0.2 GeV brick for three different times. Panel (a): 
Gluons radiated from a 16 GeV parent quark, at times 0.5 fm/c (red), 1.0 fm/c (blue), and 1.5 fm/c (green). The solid line 
shows the results of our full calculation, the dashed line the result of multiple soft scatterings, and the dotted line shows the 
results of purely elastic scatterings. Panel (b): Same as panel (a), but for a 50 GeV parent quark. Times are 0.8 fm/c (red), 
1.6 fm/c (blue), and 3.2 fm/c (green). In both panels, results obtained with AMY are shown as a black solid line. 
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FIG. 5: (Color online) Degradation of quark transverse momentum as a function of time for two QGP bricks at T — 0.2 GeV 
(panel (a)), and T = 0.4 GeV (panel (b)). Both initial distributions are delta-functions at 16 GeV and are not shown. The 
three distributions correspond to three different final times. From right to left those are: 1, 3, and 5 fm/c for panel (a), and 
0.5, 1.5, 2.5 fm/c for panel (b). 



constant we used a s = 0.3, which is similar to the value found from experiment in [7j based on the AMY framework 5 . 
Figure [3] is a repetition of this exercise with different kinematics: a 16 GeV quark now radiates a 8 GeV gluon. In 
both cases, note that the rates obtained here smoothly interpolate between those obtained at leading order in the 
opacity expansion and those produced by AMY. 

The N = 1 approximation, Eq. ((5J, is seen to accurately track the full result at small times. It begins to deviate at 
roughly the formation time, due to the onset of the Landau-Pomeranchuk-Midgal (LPM) interference effect between 
multiple scatterings occurring in the medium. The difference between the N = 1 and full curves precisely accounts 
for that effect. The amount of that difference is proportional to the number of collisions occurring during a formation 
time and can thus be used to estimate it: for RHIC energies, we thus estimate that 1 < N < 3, typically. Higher 
orders in the opacity expansion have been discussed in [221 ] ; we have not compared with these results. 



In fact, preliminary investigations suggest that the phenomenology of one-body observables - like Raa - w iU require but a modest 
increase (~15%) of a s [33j . The relative robustness of single-particle distributions is also apparent from Fig. [3] 
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The linear rise of rates at small time is thus to be interpreted within the N = 1 approximation. At small times, 
the virtuality of the jet is very high and only very hard collisions contribute. In this region, ((SJ) may be estimated as 

dT/dk » P^/p / C(q) ex P«V"/P <&in (19) 

where q 2 lin ~ fc/r and n is a number density of charge carriers. This produces the linear rise. This exhibits that the 
relevant collisions are harder than the plasma scale, and therefore more weakly coupled, making this region a robust 
prediction of QCD. 

In Fig [4] we indicate the radiation spectrum k 2 dT/dk as a function of k in a brick of temperature T = 0.2 GeV, 
for a 16 GeV parent quark with the same parameters as above. As time increases, the curves move upwards and 
eventually saturate the (time-independent) AMY prediction. In the left panel, this saturation occurs at approximately 
t = 3.5 fm/c. To avoid making the figures too busy, we have omitted the N = 1 spectra. The N = 1 spectra are 
indistinguishable from the full curve when the formation time effects are important, on the right-hand side of the 
figures (below the AMY curve), but they keep growing and overshoot the AMY asymptote on the left in keeping with 
Figures 2 and 3. 

Regarding multiple soft approximation, it is important to distinguish between the small and large time regions. 
At large times the approximation performs rather well as it tracks closely the AMY result. It produces slightly red 
spectra as may be seen in Figure[3]but also in Figure 2, where the oscillator curves asymptote to slightly higher values. 
The oscillator curves depend on a phenomenological parameter q q which we have determined, for both temperatures, 
by requiring the correct radiation rate when k — 8. That is, the asymptotic agreement in Figure [3] is a consequence 
of our choice of q q . The values we have used are q q = 0.37 GeV 2 /fm and 2.45 GeV 2 /fm at T = 0.2 GeV and T = 0.4 
GeV, respectively. These values should be compared with the definition of q q (fl~4"l) . For the collision kernel Eq. (|12l) , 
the latter evaluates to 

9 W D rio g (i + ir) 

* s — (20) 

Plugging in g max ~ 3 and 4 GeV respectively (estimated from twice the values [twice is for 2 = p/k] entering the 
transverse momentum spectra in Figure [5]) leads to q q « 0.33 GeV 2 /fm and 2.0 GeV 2 /fm respectively. We conclude 
that, in the large t region, the values of q required by the oscillator model may be larger than those corresponding to 
the quark momentum broadening problem but only by a small amount < 20%. 

In the small t region, the multiple soft scattering approximation performs poorly and does not reproduce the full 
result (which it is supposed to approximate). This is seen in Figures [5] and [3] at small times, and is also is also obvious 
from Figure S] where the ultraviolet part of the spectrum is absent at small times. In general, whenever the "Full" and 
"AMY" curves differ, meaning that interference between vacuum and medium radiation occurs, the corresponding 
suppression is exaggerated by the multiple soft scattering approximation. Analytically, at fixed p and k, the rate in 
the multiple soft scattering approximation grows like t 3 at small times while the correct, N = 1, rate rises like t. 6 
This qualitative feature has been anticipated in the above-mentioned work [2?], [28[ . Finally, we note that in all regions 
our "multiple soft approximation" curves differ by only a very small amount from those obtained within the more 
familiar massless approximation Eq. (|16p . All our conclusions thus apply also to the latter, simpler formula. 

Our choice to plot the combination k 2 dT/dk is motivated by the fact that this provides a good estimate for 
dE/dt = J kdkdT /dk, which is the relevant concept to discuss the soft region. In particular, this quantity is infrared- 
safe. At larger values k > 2 -=- 3 GeV (corresponding to the lateral shift of the measured leading-hadron spectra 
dN ~ dp±/p 8 [ 1 for Raa ~ 0.2), dE/dt ceases to be a useful concept and the more appropriate quantity is kdT/dk. 
We have included on the figure the leading-order elastic rate 

dT e \ = 9 2 m 2 D C F 

dk I6irk 2 [ ' 

as extracted from [34[ . Note that there is no kinematic cutoff on the amount of longitudinal energy which can be lost 
through elastic scattering. To the best of our understanding, elastic processes are not interference-suppressed at small 
times and so, in the brick problem, produce a time-independent contribution. As a function of the temperature, we 
find that the elastic component dominates the radiative one for k < 2.5T. 



The familiar linear rise dE/dt oc g 2 qt within the multiple soft scattering approximation |i4ll is obtained from this slow rise dT/dk oc i 3 
at fixed k combined with the rapid motion to the right of the peak value k <x qt 2 . As is evident from the figure, in RHIC kinematics 
the multiple soft scattering contribution to dE/dt is swamped out by the N = 1 contribution which rises linearly as a whole. 
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FIG. 6: (Color online) Transverse momentum distribution of 3 GeV and 8 GeV gluons radiated from a 16 GeV parent quark. 
Panel (a) is for a QGP brick at T = 0.2 GeV. From right to left the gluons are emitted at times 1, 2, 3 and 5 fm/c. Panel (b) 
is for a QGP brick at T = 0.4 GeV and times 0.5, 1 and 2 fm/c. 
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FIG. 7: (Color online) Transverse momentum distribution of a 10 GeV gluon radiated from a 50 GeV quark at T — 0.4 GeV 
and T = 0.5 GeV: Conditions of some relevance for the LHC. From right to left, times are 0.5, 1, 2, and 3 fm/c. 

This may be interpreted in two ways. First, this means that the elastic component is only important in the soft 
region, which is the standard lore. Second, this implies that soft radiation is largely irrelevant. Indeed, since there is 
no invariant way to distinguish soft emission processes from collisional ones, the best way to look at radiation with 
k -C 2.5T is as a small, sub-leading correction to the elastic component. This is relevant because (JBJ) was never meant 
to be accurate in this region anyway. This is also why we have not included medium absorption processes, as well as 
Bose-stimulation or Pauli-blocking factors, in (fT7)> - all these processes account for at most a small increase in the 
elastic component. 

In addition, it is instructive to study the efficiency with which the medium knocks apart the initial delta function, 
starting for example from a 16 GeV quark. This is shown in the two panels of Figure [5l for two different brick 
temperatures. At the later times, any hint of a peak structure has mostly been washed out. Notably, the inclusion 
of the formation time effects does reduce the distribution in magnitude, but does not introduce any significant shape 
distortion. 7 

To address the self-consistency of the formalism it will be important to look at the k± distribution of the radiated 
gluons. This is shown in Figure |fJ] It is difficult to draw any firm conclusion regarding this issue in the RHIC context 
at present since in this paper we have focused on a semi-realistic brick problem. Nonetheless, in the cases indicated we 





T=0.4 GeV — 






T=0.5 GeV 













7 In producing the AMY curves in Figure [5] radiation was included only for t > 0.7 fm. Comparable cutoffs were used in previous AMY 
implementations ;7j, where they coincide with the cutoff on the onset of hydrodynamics, and the resemblance of this cutoff with the 
formation times in Figures 2 and 3 largely explains the apparent smallness of the effect. 
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are satisfied that the relevant values of k± are always smaller than the energy k of the gluon. In the worst case we have 
considered, corresponding to a 3 GeV gluon radiated from a 16 GeV parent quark in a T = 0.4 GeV brick, the typical 
opening angle from k±/k is approximately 30°. In other cases it is much smaller. The collinear approximation should 
thus be under control in the cases shown in this paper. The k± scale is mostly determined by the temperature of the 
medium, although k also plays a role (this is in qualitative agreement with the parametric estimate k± ~ (gfc) 1 / 4 in 
the deep LPM regime). The kj_ scale in the figure gives the scale at which the running coupling is to be evaluated 
(up to an unknown multiplicative factor whose determination would require a full 1-loop computation), which should 
be useful for the reliable inclusion of running coupling effects in the future. In many cases the fcj_'s relevant at RHIC 
are seen not to be much larger than 1 GeV, suggesting a possible importance of loop corrections. On the other hand, 
strong non-perturbative effects would appear unlikely at such transverse momentum scales. The situation on all these 
issues is expected to improve under LHC conditions, as shown in Figure [7] 

VI. CONCLUSION 

We have described a method to solve exactly the BDMPS-Z radiation formula. The key step in this method is the 
focus on a differential radiation rate Eq. ((6]), which is both numerically tractable and physically meaningful. Such 
a rate formulation immediately connects with the AMY formalism. On the other hand, this rate is constructed in 
such a way that its time integral agrees precisely with the radiation probability previously defined by BDMPS-Z, 
including in situations in which the duration of the emission process is not negligible and where the AMY treatment 
was previously not applicable. This includes interference effects between medium and vacuum radiation. The method 
makes it possible to study numerically and exactly, for the first time, BDMPS-Z radiation spectra in all kinematic 
regimes. The rate formulation makes the transitions between different regimes particularly easy to see. 

In the brick problem, we have distinguished two different interference-dominated regimes. In the first one, occurring 
at small time, the medium modifications to the radiation are determined by interference between the vacuum radiation 
and that produced by a single hard collision. This regime smoothly connects with a second one, in which the vacuum 
radiation becomes irrelevant and the radiation is determined by Landau-Pomeranchuk-Migdal (LPM) interference 
between in-medium radiation. The latter regime is well described both by AMY and by the multiple soft scattering 
approximation, while the former regime is well described by the GLV N = 1 approximation. The transition between 
these two regimes occurs at a physically important length scale, t ~ 1 -j- 2 fm depending upon the energies, which 
makes it essential to accurately describe it. 

Generally speaking, for a given medium we have shown that AMY (in previous implementations) over-estimates 
radiation while the multiple soft scattering approximation systematically under-estimates it. This agrees qualitatively 
with the pattern of interaction rates which have been extracted from experiment in the past using those models [7] . 

Phenomenologically speaking, information about the underlying medium enters exactly two places in the present 
formalism. One is the elastic collision rate d 2 T c \(q)/d 2 q, related to the quantity C(q) in this paper, and the other is 
the thermal masses. Although we have made no detailed investigation of the effect of thermal masses it seems clear 
that the most important parameter is the elastic collision rate. For a plasma isotropic in its rest frame, this is a 
function of |q| which is independent of the jet energy. Our analysis shows that its second moment, q q , does not suffice 
to characterize jet energy loss as it ignores a hard Coulomb tail which dominates at small times. In the present paper 
we have contented ourselves with the simplest available choice, Eq. (|12[) , but it will be important in the future to 
learn more about the phenomenological implications of the shape of this curve. 

It is known that current efforts to extract meaningful physics from jet quenching measurements in heavy ion collisions 
have to be done with the help of hydrodynamic modelling. The formation time effects discussed in the present paper 
are expected to be particularly relevant for non-central collisions, and for angular correlation measurements, where 
the details of the geometry of the emitting region play a larger role. As a consequence, a systematic analysis of several 
different observables will be needed to quantitatively characterize the hot and dense strongly interacting medium 
formed at RHIC, and soon to be created at the LHC. 
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Appendix A: Computing the radiation in an expanding medium 

In this appendix, we spell out, for future reference, the proper formulae to be used in expanding media. The 
formula will also be valid in the situations in which the medium is not homogeneous in the longitudinal direction over 
a formation time. The derivation follows the discussion in section IIII Al where longitudinal homogeneity was indeed 
never required. 

The generalization of ([5]) in this case is 

^ = ^M xRe /V/ =JS—Mt)C(t)K(t,c l ;t 1 ,p).p. (Al) 
dk np J Jq lP SE(q,t) 

In this expression i/5E(q,t) = ffdt'e-'fi SE . (This will always be very close to i/SE(q, t), provided thermal masses 
change slowly.) 

Equation © contains the time-dilation factor ^(t) = g ^ a3ma , It describes the transformation of d 2 r c i/d 2 q from the 
plasma frame to the lab frame. This factor always multiply the operator C. If, in the lab frame, the plasma is moving 
with velocity v and at angle 9 relative to the jet axis, this factor may be written 

1 — v cos 6 

The expression Q for the now time-dependent C is not otherwise modified and C (q) entering it should be computed 
in the plasma frame as usual. 

The same factor multiply C in the equation for K which now reads 

[d t + iSE(q, t) + 7(*)C(t)] K(t, q; h, p) = 0, K(h , q; h, p) = -J- , (2n) 2 6 2 (q - p). (A2) 



Ak 2 {p-k) 2 



The expression for the light-cone energy is unchanged 



*FVn t\ - I | ™M m a^) (AO's 

^ (q '* } - 2k{p-k) ~W 2^fc) ~ ^2p~ (A3) 

provided all energies in it are measured in the lab frame. For heavy quarks, the masses are the rest masses, while for 

light quarks and gluons the masses are the asymptotic thermal masses (respectively = 3 and = — in 
QCD at leading order in thermal perturbation theory). 

The solution to these equations exhibits, in an evolving medium, in addition to the finite-size effects discussed in 
the present paper, important memory effects. Their analysis lies outside our present scope. 

We briefly describe our numerical implementation of (lAip . We view K(t\ t{) as an evolution operator in t\ starting 
from time t. (Since t\ < t, the evolution time is \t\— The initial wavefunction at time t is the vector-valued 
function ^(p, |ti— £|=0) = j(t)C(t) o =2—. We evolve it numerically in the "interaction picture" as 

MP, \ti-t\) = -e*& SE <» ) m)C(t 1 )e- i & SE{P) MP, \h-t\) 



d\h-t\ 



Going to the interaction picture removes violent phases from ipj. The initial conditions at time t decay at large 
p like ipi ~ pC(p) ~ p/p 4 - This decay is preserved at all times. To find T we need to take the integral 
J J dtye 5E ( p ^ p-^/(p, \t\— 1\). Because the p integral is not absolutely convergent, it is most convenient to 
perform the t\ integration first. At large times, tpi tends to zero and the time integral saturates. However, oscillations 
in time will cause an extra suppression ~ l/SE(p) at large p making the p integration convergent, as it should be 
expected physically. Rotational invariance in the p plane can be used to set tpi = -T-?/>/(|p|) for some scalar function 
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